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Abstract 

In this paper we continue our development of a dimensional perturbation theory (DPT) treatment of A'^ 
identical particles under quantum confinement. DPT is a beyond- mean-field method which is applicable 
to both weakly and strongly-interacting systems and can be used to connect both limits. In a previous 
paper we developed the formalism for low-order energies and excitation frequencies. This formalism 
has been applied to atoms, Bose-Einstein condensates and quantum dots. One major advantage of the 
method is that A'' appears as a parameter in the analytical expressions for the energy and so results for N 
up to a few thousand are easy to obtain. Other properties however, are also of interest, for example the 
density profile in the case of a BEC, and larger A'' results are desirable as well. The latter case requires 
us to go to higher orders in DPT. These calculations require as input zeroth-order wave functions and 
this paper, along with a subsequent paper, address this issue. 



1 Introduction 

Systems that involve N identical interacting particles under quantum confinement appear throughout many 
areas of physics, including chemical, condensed matter, and atomic physics[l]. In a previous paper[2], we 
applied the methods of dimensional perturbation theory (DPT) [3] , a powerful set of tools that uses symmetry 
to yield simple results, to such iV-body systems. We presented a detailed discussion of the dimensional 
continuation of the iV-particle Schrodinger equation, the D ^ oo equilibrium (-D'^) structure, and the (D^^) 
energy. We used the FG matrix method to derive general, analytical expressions for the A^-body normal- 
mode vibrational frequencies, and we gave specific analytical results for three confined A^-body quantum 
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systems: the AT-electron atom, the A''-electron quantum dot, and the A''-atom inhomogeneous Bose-Einstein 
condensate with a repulsive hardcore potential. 

In a subsequent paper [4] we further pursued our analysis of the ./V-atom inhomogeneous Bose-Einstein 
condensate (BEC)with a repulsive hardcore potential and optimized our low-order DPT energy by fitting to 
low-N DMC data[5] and then extrapolated the result out to large N. Our low-order energies are accurate 
out to -/V w 10, 000 and larger N could in principle be obtained by going to higher orders in the perturbation 
expansion [6]. Ten thousand atoms may not seem particularly impressive given the fact that the Gross- 
Pitaevskii (GP) equation can yield higher-A'' results, but the advantage of DPT is that it is not a weak- 
interatomic-interaction approximation. At lowest orders it includes beyond-mean-field effects. Thus DPT is 
effective for large-a systems that have been created in the lab by exploiting Feshbach resonances [7], systems 
for which the GP and modified Gross-Pitaevskii (MGP) equations fail. This we demonstrated for a BEC 
which had an a one hundred times that of the natural scattering length of ^''Rb. 

In another paper [8], low-order A/'-body dimensional perturbation methods have been applied to the A''- 
electron atom. In this instrumental paper, Loeser obtained low-order, analytical expressions for the ground- 
state energy of neutral atoms. For Z = 1 to 127, the numerical results compare well to Hartree-Fock energies 
with a correlation correction. 

All of the previous work has focussed on low-order calculations of energies and thus far little attention 
has been given to wave functions, even at lowest order. Yet the lowest-order wave function yields important 
information about the system. The normal mode coordinates tell us the nature of the excitations of the 
system and the lowest-order wave function gives us expectation values at low orders. For macroscopic 
quantum-confined systems, such as the BEC, the wave function is uncloaked in an explicit way since the 
density profile may be viewed in a direct fashion. Also, calculating energies and wave functions to higher 
orders in 1/Z) requires as input the lowest-order wave function. 

We briefly outlined the derivation of the lowest-order wave function in a previous letter [9] and in this 
paper we expand considerably upon that discussion. In Section 2 we discuss the S-wave, Jacobian-weighted 
Hamiltonian and wave function. In Section 3 we dimensionally scale the system to regularize the large 
dimension limit for which we then solve. In Section 4 we begin the normal mode analysis which will result 
in the energy through first order in 1/D and the lowest-order wave function. In Section 5 we discuss 
the symmetry of the A''-body quantum-confinement problem. This symmetry simplifies the problem to an 
extraordinary degree, making an exact solution of this ./V-body problem, with N{N — l)/2 interparticle 
interactions, possible at large D . In Section 6 we exploit this symmetry further and discuss the use of 
symmetry coordinates in the calculation of the large- dimension, zeroth-order normal modes of the system. 
In Section 7 we introduce a particular approach to determining a suitable set of symmetry coordinates. 
The method uses what we term primitive irreducible coordinates, from which the symmetry coordinates are 
derived. This is illustrated in detail in a simple example involving the coordinates of a particular species of 
the A''-body quantum-confinement problem (symmetry coordinates transforming under the same irreducible 
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representation are said to belong to the same species). In Section 8 we apply the general theory discussed 
in earlier sections to derive in detail the frequencies and normal-mode coordinates of this species. Section 9 
is a summary and conclusion section. The frequencies and normal-mode coordinates of the remaining, and 
more complicated species of the system are derived in a subsequent paper [10]. 



2 The D-dimensional N-body Schrodinger Equation 



For an A''-body system of particles confined by a spherically symmetric potential and interacting via a 
common two-body potential gij, the Schrodinger equation in D-dimensional Cartesian coordinates is 



N N-1 N 

i=\ i=l j=i-\-l 



where 



D 



2m, ^ dxl 



conf 



and 



9i3 = Knt ^^Xl^^i {Xi„-Xj^f^ 



(1) 



(2) 



(3) 



are the single-particle Hamiltonian and the two-body interaction potential, respectively. The operator H is 
the Z)-dimensional Hamiltonian, and is the u*^ Cartesian component of the i*'* particle. The term Vconf is 
the confining potential. For the A''-electron atom, Fconf is the Coulomb attraction between the electrons and 
the nucleus, while for the ^''-electron quantum dot and ./V-atom hard-sphere problem, a model appropriate 
for trapped gaseous BECs, we model the confinement as a harmonic trapping potential. The two-body 
interaction potential Vint is Coulombic in the first two systems and a hard sphere in the third. 



2.1 Transformation of the Laplacian 



Restricting our attention to spherically symmetric [L = 0) states of the iV-body system, we transform from 
Cartesian to internal coordinates. A convenient internal coordinate system for confined systems is 



D 



{l<i<N) and 7^- 



j) = Y^Xi^Xj^ /rirj {l<i<j< N), 



(4) 



which are the D-dimensional scalar radii r, of the N particles from the center of the confining potential and 
the cosines 7,j of the N{N — l)/2 angles between the radial vectors. 

Now for a function ^ dependent on only two functions r{x) and 7(3;) one can write 



(P'S{r{x),j{x)) d^r ^^7 f dr\^ d'^'i 



dx^ 



dx'^ dr dx"^ d'y \rfa; 



dr^ 



+ 



2 J2 



(i^* ^d'y dr d^* 
l_ 2 - . 

^7^ dx dx drd-y ' 



(5) 
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Generalizing this, when operating on the state ^{ri{xii,),'yii{xi„) . . .^ikixiv) ■ ■ ■JiNixiv)) (where k ^ i and 
v = 1,. . . ,D)), Vf can be written in terms of the internal coordinates of Eq. (4) as 

\dxiv J \dxiu) dvid^ij 
The relevant derivatives of the internal coordinates are 



which lead to the effective 5-wave Laplacian in internal coordinates: 

D-1 d ^ x-^D-lx-^ d 



^7.. (7) 



2.2 The Jacobian- Weighted Hamiltonian 

All of the -/V-body DPT work up till now has utilized a similarity transformation so that the kinetic energy 
operator is transformed to a sum of terms of two kinds, namely, derivative terms and a repulsive centrifugal- 
like term, which when attractive interparticle potentials are present, stabilizes the system against collapse 
in the large-Z) limit. The particular transformation which has been used in almost all the work to date, is 
chosen so that first-order derivatives resulting from the Laplacian (Eq. (9)) are removed from the transformed 
Hamiltonian. When this is done the zeroth and first orders of the dimensional (1/-D) expansion of the 
Hamiltonian become exactly soluble for any value of A''. In the — > oo limit, the second derivative terms 
drop out, resulting in a static problem at zeroth order, while first order corrections correspond to simple 
harmonic normal-mode oscillations about the infinite-dimensional structure. 

Writing the similarity transformation of the wave function ^ and operators O as 

$ = X"'*, and O = x'^Ox, (10) 

the transforming function is chosen to be 

X=(r-ir2...r^)-(^-^)/^r-(^-i)/^ (11) 

where F is the Gramian determinant, the determinant of the matrix whose elements are 7ij (see Appendix B). 
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Carrying out the transformation of the Hamiltonian and wave function of Eq. (1) via Eqs. (10) and (11), 
the Schrodinger equation becomes [8]: 



((i»)T^+l^)(^^^)* = -E(,ii)$ (12) 




Af-l N 

V = J2 ^conf (rO + E E V^r^tinj). (14) 

i=l i=l j=i-\-l 

The Gramian matrix whose determinant is F^*^ is the i*^ principal minor formed 1)\ deleting from F the row 
and column corresponding to the i*^ particle. The quantity nj = y^r? + r| — ^riVj^ij is the interparticle 
separation. The similarity transformed Hamiltonian for the energy eigenstate (iu)^ is (m)!! = x~^iJx, 
where (m)!! = {^^^^■^T + V). As we can see from Eq. (13) all first-order derivatives have been eliminated from 
the Hamiltonian. 

While the similarity transformation of Eqs. (10) and (11) is convenient in that all first-order derivative 
terms are eliminated, it complicates the interpretation of the normal-mode structure of the large-!) wave 
function. Only when the weight function, W, for the matrix elements is equal to unity, i.e. 

W = Jx^ = l (15) 

or 

X = J-"\ (16) 

where J is the Jacobian of the transformation to internal coordinates, will a first derivative of an internal 
coordinate be the conjugate momentum to that coordinate. The transforming function of Eq. (11) does not 
satisfy Eq. (16) since the Jacobian J is[ll] 

J=(rir2...r^)^-iF(^-^-i)/2. (17) 

Furthermore, the development of higher-order DPT expansions involves the matrix elements of coordinates 
and their derivatives between the zeroth-order normal-mode functions. These matrix elements are much 
easier to calculate when the weight function in the integrals is unity, and this only occurs when Eq (16) is 
satisfied. 

In Ref. [11], Avery ct al. considered the problem of performing the similarity transformation of Eq. (10) 
with a more general transforming function which has adjustable parameters a and /3, and is of the form: 

X=(rir2...rAr)-"F-^/'. (18) 

One of the cases considered by Avery et ai.[ll] has a = {D — l)/2 and (5 = {D — N — 1) /2 so that 

X = {nr2 . . . r-;v)-^^-'^/'F-(^-^-i)/4 . (19) 
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The transforming function x of Eq. (19) does satisiy Eqs. (15) and (16). 

Carrying out the transformation of the Hamiltonian and wave function of Eq. (1) via Eqs. (10) and (19), 
the Schrodinger equation becomes[ll]: 



where 



N 



T = h^y: 



i=l 



1 ^2 



+ 



2mi dn^ 
Ar(Ar - 2) + (£> - TV _ 1)2 (^rw ^ 



- o ^ 2 ( y y^bjk - lijlik)-^-^ NY] 



911 



N 



Snurf 



1 1 v-v- 5 . ^ 9 

>2 fr'') 



2mi dri^ 2mirf 



ik 



N{N-2) + {D-N-iy 



+- 



(20) 



(21) 



The latter expression for (j-jT is explicitly self-adjoint since the weight function, W, for the matrix elements 
is equal to unity, i.e. Eq. (15) is satisfied. The similarity transformed Hamiltonian for the energy eigenstate 

(i)$is (,)ff = ((,)r + y). 



3 Infinite-D analysis: Leading order energy 

To begin the perturbation analysis we regularize the large-dimension limit of the Schrodinger equation by 
defining dimensionally scaled variables: 

fi = ri/K{D) , E = k{D)E and H = k{D) (^^H (22) 

with dimension-dependent scale factor k{D). Prom Eq. (21) one can see that the kinetic energy T scales in 
the same way as 1/r^, so the scaled version of Eq. (20) becomes 

= {-^f + Kffl * = E^, (23) 



AD)' 

where barred quantities simply indicate that the variables are now in scaled units. Because of the quadratic 
D dependence in the centrifugal-like term in T of Eq. (21), we conclude that the scale factor k{D) must also 
be quadratic in D, otherwise the Z) — > oo limit of the Hamiltonian would not be finite. The precise form of 
k{D) is chosen so that the result of the scaling is as simple as possible and depends on the system in question. 
In previous work[2] we have chosen k{D) = {D — 1){D — 2N — 1)/(4Z) for the 5- wave, A''-electron atom; 
k{D) = QZho for the AT-electron quantum dot where Q = {D — 1){D — 2N — l)/4 and the dimensionally- 
scaled harmonic oscillator length and trap frequency respectively are Zho = y ra'an, '^^o = ^^'^ho > ^"^^ 



k{D) = D'^a^a for the EEC where a^o = \ '^o = D^i^o ■ The factor of k{D) in the denominator of 
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Eq. (23) acts as an effective mass that increases with D, causing the derivative terms to become suppressed 
while leaving behind a centrifugal-like term in an effective potential, 



N 



Veff(r,7;<5 = 0)=^ 



8mif2 r 



N-1 N 



+ v;onf(f,7;<5 = o) +^ ^ v'int(f,7;<5 = o), 



(24) 



i=l ^ ■- ^ ' i=l j=i+l 

where 5 = l/D, in which the particles become frozen at large D. In the — > oo (5^0) limit, the excited 
states have collapsed onto the ground state, which is found at the minimum of V^tt . 

We assume a totally symmetric minimum characterized by the equality of all radii and angle cosines of 
the particles when Z) — > oo, i.e. 



Ti = r„ 



{l<i<N), iij=l^ {l<i<3<N). 



(25) 



Since each particle radius and angle cosine is equivalent, we can take derivatives with respect to an arbitrary 
fj and 7ij in the minimization procedure. Then evaluating all fj and 7,j at the infinite-f radius and angle 
cosine, foo and , respectively, we find that foo and satisfy 



dfi 

dVeff (r,7;(5) 



(5=0 



5=0 



= 



= 0, 



(26) 
(27) 



where the 6 = subscript means to evaluate all r, at Voo and all jij at 'y^. In scaled units the zeroth-order 
{D ^ oo) approximation for the energy becomes 



Eoo = V;ff(rcx,,7oo;0)- 



(28) 



In this leading order approximation, the centrifugal-like term that appears in Vetf , even for the ground state, 
is a zero-point energy contribution required by the minimum uncertainty principle [12]. 



4 Normal- mode analysis and the 1/D first-order quantum energy 
correction 

At zeroth-order, the particles can be viewed as frozen in a completely symmetric, high-£) configuration 
or simplex, which is somewhat analogous to the Lewis structure in atomic physics terminology. Likewise, 
the first-order 1/D correction can be viewed as small oscillations of this structure, analogous to Langmuir 
oscillations. Solving Eqs. (26) and (27) for foo and gives the infinite-!? structure and zeroth-order energy 
and provides the starting point for the 1/D expansion. To obtain the 1/D quantum correction to the energy 
for large but finite values of D, we expand about the minimum of the D — > oo effective potential. We first 
define a position vector, consisting of all N{N + l)/2 internal coordinates: 
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where 



and 



r-2 



(30) 



7 




7l5 
725 



(31) 



\ IN-I.N ) 

We then make the foUowing substitutions for all radii and angle cosines: 

f,=foo+5^/M (32) 

7ii=7oo+'^'/'%-, (33) 

where 5 = 1/2? is the expansion parameter. Now in all practical situations Kff (^; S) is a function of y and 
5 and so we may then obtain a power series in S^l'^ of the effective potential about the Z) — > oo symmetric 
minimum as 



Veff(y;5) = 









■ d8 ' 


dS 












51/2=0 





51/2=0 



E 













. dyt. 


51/2. 


51/2=0 


a<5i/2 


51/2=0} 



4- 



9Veff 



aVeff 






■ d^s ■ 


+ 








■ dS ' 


5(5 




51/2=0 


(951/2)2 


51/2=0 


5(52 


I'm. 


51/2=0 





+ 



51/2=0 



dy^L 
p p 

+EE 

I^— 1 



_5 
55 



51/2 



dyt. 







' dy^ ' 




■ 5,5 ■ 






51/2=0 


.5(51/2 


51/2=0 


5(51/2 


51/2=0 



51/2=0 



(5,51/2)2 



51/2=0 







[ 52v;if 1 




' 5y^ ■ 




_5,5i/2_ 


51/2=0 


dyfj^dy^ 


51/2=0 


5,51/2 


51/2=0} 



+0(53/^), 



(34) 



where 



P = iV(iV + l)/2 



(35) 
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is the number of internal coordinates. The 0((^^/^)'^) term in the power series (Eq. 43) is simply the 
zeroth-order energy (Eq. 28). The 0((^^/^)^) term is zero since 

dS 



= 



(36) 



51/2=0 

and we are expanding about the minimum of the effective potential so 

= 0, 











51/2. 


51/2=0 



(37) 

i.e. the system is said to be in equilibrium since the forces acting on the system vanish [Eqs. (26) and (27)]. 

Defining a displacement vector consisting of the internal displacement coordinates [primed in Eqs. (32) 
and (33)] 

^ f 
7' 



where 



and 



we then have 



and also defining 



Eq. (34) becomes 



y 



7 




V 7iv-i,jv / 



dyf. 



51/2=0 



85 



Vi^A 51/2=0 



P P 



v;« (y'; ^) = ,v2=o + J M E E 



2 l^y^^^y-^ 



(38) 



(39) 



(40) 



+ Vo> + 



o 



51/2=0 



(41) 



(42) 



(43) 
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The first term of the 0(((5^/^)^) term defines the elements of the Hessian matrix[13] F of Eq. (45) below. 
The derivative terms in the kinetic energy are taken into account by a similar series expansion, beginning 
with a first-order term that is bilinear in d/dy' , i.e. 

^ = -\^jijiG,,dy,^dy. + O (5^'^) , (44) 

where T is the derivative portion of the kinetic energy T (see Eq. (21)). Thus, obtaining the first-order 
energy correction is reduced to a harmonic problem, which is solved by obtaining the normal modes of the 
system. 

We use the Wilson FG matrix method[14] to obtain the normal-mode vibrations and, thereby, the first- 
order energy correction. It follows from Eqs. (43) and (44) that G and F, both constant matrices, are 
defined in the first-order S = 1/D Hamiltonian as follows: 

Hi = - \dy'^Gdy, + \y'^Fv' + Vo . (45) 

After the Schrodingcr equation (20) has been dimensionally scaled, the second-order derivative terms are of 
order 5, and by comparing these terms with the first part of Hi , the elements of the kinetic-energy matrix 
G are easily determined. The elements of the Hessian matrix[13], F^j, = ^ ^ _^ , on the other hand, 
require a bit more eflfort to obtain [2]. 

We include a derivation of the FG matrix method in Appendix A, but we state here the main results of 
the method. The h^^ normal mode coordinate may be written as (Eq. (A19)) 

[q']b = h^y\ (46) 
where the coefiicient vector h satisfies the eigenvalue equation (Eq. (A20)) 

VGb = \b (47) 

with the resultant secular equation (Eq. (A21)) 

det(FG - AI) = 0. (48) 
The coefiicient vector also satisfies the normalization condition (Eq. (A22)) 

b^Gb = l. (49) 
As can be seen from Eq. (A5) in Appendix A the frequencies are given by 

Af, = qI, (50) 
while the wave function is a product of P = N{N + 1) /2 harmonic oscillator wave functions 

^o{y') = f[hn,{u>l/'[q%) , (51) 
6=1 
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where (j^l^^Wlb^ is a one-dimensional harmonic-oscillator wave function of frequency Qb , and nt is the 
oscillator quantum number, < < oo, which counts the number of quanta in each normal mode. 

Depending on the number of particles, the number of roots A of Eq. (48) - there are P = N{N + l)/2 
roots - is potentially huge. However, due to the ^jv symmetry (see Ref. [15] and App. C) of the problem 
(see Sect. 5), there is a remarkable simplification. Equation (48) has only five distinct roots, A^, where fi 
runs over 0~, 0+, 1~, 1+, and 2, regardless of the number of particles in the system (see Refs. [2, 8] and 
Sect. 6.3). Thus the energy through first-order (see Eq. (53)) can be written in terms of the five distinct 
normal-mode vibrational frequencies, which are related to the roots A^ of FG by 



A^ = . 



The energy through first-order in S = 1/D is then [2] 



E = E^+5 



n„=0 



(52) 



(53) 



2 } 



where the are the vibrational quantum numbers of the normal modes of the same frequency tZi^ (as such, 
counts the number of nodes in a given normal mode). The quantity d^,„^ is the occupancy of the manifold 
of normal modes with vibrational quantum number and normal mode frequency ui^, i.e. it is the number 
of normal modes with the same frequency and the same number of quanta n^. The total occupancy of 
the normal modes with frequency u)^ is equal to the multiplicity of the root A^, i.e. 



— E ^/^'"i* 



(54) 



where is the multiplicity of the /x*" root. The multiplicities of the five roots are[2] 



do+ 
do- 
di+ 

d2 



1, 
1, 

N-1, 
N-1, 
N{N-3)/2. 



(55) 



Note that although the equation in Ref. [8] for the energy through Langmuir order is the same as Eq. (53), 
it is expressed a little differently (See Ref. [16]). 



5 The Symmetry of the Large-D, N-body Quantum-Confinement 
Problem 

Such a high degree of degeneracy of the large-D, A^-body quantum-confinement problem indicates a high 
degree of symmetry. In this section we study this in more detail. 
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5.1 The Indical structure of F, G, and FG matrices 



The F, G, and FG matrices, which wc gcucrically denote by Q, are P x P matrices with the same indical 
structure as yy"^: 





r'lr', ■ 


• r[f'N 






r'll'N-iN 




r'2r'2 ■ 




r'21'12 


^2713 


r'21'N-lN 






^'n^'n 


r'Nl'i2 




r'nl'N-iN 


7'l2^'l 


7'l2^2 • 


I'ur'N 


7'l27'l2 


7'l27i3 


1'i21'n-\n 




713^2 • 


7'l2^1 


7i37i2 


7'l37'l3 


I'izI'n-in 








7jv-ijv7'i2 


7jv-ijv7'i3 • 


i'n-in1'n-in 



(56) 



where y is defined by Eq. (29). The indical structure of this matrix suggests a convenient shorthand for 
referencing the elements of the Q matrices. The upper left block of Eq. (56) is an (A'' x A'') matrix with 
elements associated with {f\,f'j)\ hence we use the subscript (i, j) to refer to these elements. The upper right 
block is an {N x N{N — l)/2) matrix with elements associated with (f-,7j-;.); hence, we use the subscript 
{i,jk) to refer to these elements. The lower left block is an {N{N— l)/2xN) matrix with elements associated 
with {%j,f'i^); hence, we use the subscript {ij, k) to refer to these elements. Finally, the lower right block is 
an {N{N — l)/2 x A''(A'' — l)/2) matrix with elements associated with (7ij,7/j;); hence, we use the subscript 
{ij, kl) to refer to these elements. 



5.2 The Sn symmetry of the Q matrices 

As the number of particles N increases, diagonalizing the P x P FG matrix (where P = N{N + l)/2) 
becomes, prima facie, a daunting task. However, one of the advantages of dimensional perturbation theory is 
the simplifications that occur in the large- dimension limit. In particular, since we are dealing with identical 
particles in a totally symmetric configuration (the Lewis structure) in which all the particles are equivalent, 
the Q matrices display a high degree of symmetry with many identical elements. Specifically, 





= Qi',i' 




Qa 






Qi,j 


— Qi' ,j' 




Qb 


{^ 


^ j) and (^' ^ /) 


Qij,i 


— Qi'j'^i' 




Qc 


i^ 


^j) and (I'^f) 


Qjk,i 


— Qj'k'^i' 




Qd 


(i 


j^jj^k) and {i' ^ f ^ k') 




— Qi'^i'j' 




Qe 


(i 


^3) and ii'^j') 




— Qi'.j'k' 




Qf 


(i 


^j^k) and {i' ^ j' ^ k') 


Qij,ij 


— Qi'j',i'j' 




Qa 


i^ 


^j) and ii'^f) 


Qij,jk 


— Qi'j',j'k' 




Qh 


{^ 


^ jV fc) and (z' ^ f ^ k') 


Qij,kl 


= Qi'j\k'l' 




Q. 


i^ 


^j^k^l) and {i' ^ f ^k'^ V) 
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Note the indices in the relationships above run over all particles (1, 2, . . . , A'') with the exceptions noted in 
the far right column. For example, Qij = Qi'j' = Qb, where {i ^ j) and {i' ^ j'), means that all off-diagonal 
elements of the upper left block (the pure radial block) of Q are equal to the same constant Qb- Similarly, 
Qij,ki = Qi'j',k'i' = Qi, where {i ^ j ^ k ^ I) and {i' ^ j' ^ k' ^ I'), means that any elements of Q in 
the lower right block (the pure angular block) that do not have a repeated index are all equal to the same 
constant Q^. 

Equations (57) imply that the Q matrices, F, G, and FG, are invariant under interchange of the particle 
labels. Thus Eqs. (57) show that the system is invariant under the group of all the permutations of N objects, 
where the objects are the particle label indices. This is the symmetric group Sn (see Appendix C). As we 
shall see, this Sn invariance greatly simplifies the problem which would otherwise be intractable for large 
N; allowing us to solve for the N{N + l)/2 normal-mode coordinates, and hence the Langmuir-order wave 
function. The Sn symmetry is also behind the remarkable reduction in the number of distinct frequencies 
we noted at the end of Sec. 4, from a possible A''(A'' -|- 1)/2 to only five distinct frequencies. 

Parenthetically, we should remark here that G and F are also symmetric matrices (G-^ = G and F-^ = F); 
however, while FG does display the high degree of symmetry of Eqs. (57), it is not a symmetric matrix. 

5.3 Q matrices in terms of simple submatrices 

The Sn symmetry of the Q matrices {F, G, and FG) described in Eqs. (57) allows us to write these matrices 
in terms of six simple submatrices. We first define the number of 7- • coordinates to be 



and let Ijv be an A?^ X A/^ identity matrix, Im an M x M identity matrix, 3n 3i.n N x N matrix of ones and 
3m an M X M matrix of ones. Further, we let R be an A'' x M matrix such that Rijk = ^ij + ^ik-, 3nm be 
an A'' X M matrix of ones, and 3%^ — ^mn- 
We then write the Q matrices as 



where the block Qfff has dimension (NxN), block Qf^' has dimension (NxM), block Q^'f has dimension 
{M X N), and block Q^'^> has dimension (M x M). Now, as we show in Appendix B of Ref. [2], Eqs. (57) 
allows us to write the following: 



M = N{N -l)/2, 



(58) 




(59) 



{Qa - Qb)lN + QbJw 



(60) 



Qf'y' = {Qe — Qf)^+ Qf-^NM 
Qy'f' = (Qc — Qd}^^ + QdJ^M 

Q^.^. = {Qg-2Qh + Q,)lM + {Qh-Q>.)'R^'R + Q,3M. 



(63) 



(61) 



(62) 
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In particular, letting Q = FG, the matrix that must be diagonalized, Eq. (59) becomes 
where we have used the following abbreviations: 



d 


= 




- {FG)b = {F, - Fb)Ga 




b 






= FbGa 




c 




(FG), 


- {FG)d = (Fe - Fj)Ga 




d 




(FG)d 


= FfGa 




e 




{FG), 


- {FG)f ^ (Fe - Ff){Gg + {N- A)Gh) 




f 




iFG)f 


= 2F^Gh + Ff{Gg + 2{N - 3)G^) 




9 




(FG), 


- 2{FG)h + {FG), = [Fg - 2Fh + F,){Gg ■ 


-2Gh) 


h 




{FG)h 


- (FG), = FgGh + Fh{Gg + {N- Q)Gh) - 


- F,{Gg + [N 


I 




{FG\ 


= AFhGh + F,{Gg + 2{N - A)Gh) • 





(65) 



where 



The right-hand sides of Eq. (65), the FG matrices expressed in terms of the F and G matrix elements, may 
be derived using the graph-theoretic techniques discussed in Appendix B of Rcf. [2] . 

Wc also require the G matrix for the normalization condition (Eq. (49)). It has a simpler structure than 
the FG matrix, 

, d'lN 

G = ~ h (66) 

-g'lM + /I'R^R 

a' = {G)a 

~g' ^ {G)g-2{G)H (67) 
h' = {G)h 

and 

{G)b = 
(G)c = 

{G)d = (68) 
{G)e = 
{G)f = 
(G). = 0. 

The quantities (G)a, {G)g and {G)h depend on the choice of k{D) (See Ref. [2]). 
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6 Symmetry and Normal Coordinates 



The FG matrix is a N{N + l)/2 x N{N + l)/2 dimensional matrix (there being N{N + l)/2 internal 
coordinates), and so Eqs. (47) and (48) could have up to N{N + l)/2 distinct frequencies. However, as noted 
above, there are only five distinct frequencies. The Sn symmetry is responsible for the remarkable reduction 
from N{N + l)/2 possible distinct frequencies to five actual distinct frequencies. As we shall also see, the 
Sn symmetry greatly simplifies the determination of the normal coordinates and hence the solution of the 
large-D problem. 

6.1 Symmetrized coordinates 

To understand why there are only five distinct frequencies and to start the process of solving for the normal 
coordinates, we first look at the Sn transformation properties of the internal coordinates. 

We first note that the internal-coordinate displacement vectors f' and 7' of Eqs. (39) and (40) are basis 
functions which transform under matrix representations of Sn and each span the corresponding carrier 
spaces. These representations of Sm however, are not irreducible representations of Sm (Appendix C). 

We will show that the reducible representation under which f' transforms is reducible to one 1-dimensional 
irreducible representation labelled by the partition [N] (the partition denotes a corresponding Young diagram 
( = Young pattern = Young shape) of an irreducible representation (see Appendix C)) and one {N — 1)- 
dimensional irreducible representation labelled by the partition [A^ — 1, 1]. We will also show that the 
reducible representation under which 7' transforms is reducible to one 1-dimensional irreducible representa- 
tion labelled by the partition [N], one (A^— l)-dimensional irreducible representation labelled by the partition 
[A'" — 1, 1] and one N[N — 3)/2-dimensional irreducible representation labelled by the partition [N — 2, 2]. 

6.1.1 The Characters of the Reducible Representations of Sn under which f' and 7' Trans- 
form. 

Let us denote the cycle structure of a permutation (see Appendix C) by the symbol (l''^ , 2''^ , S''^' , . . . , N"'^), 
where the notation j"^^ means a cycle of length j and vj equals the number of cycles of length j in that 
permutation[15]. Crucially the characters, t^(i''^'2''2,...,jv''«)^ elements of a matrix representation of 

Sn corresponding to permutations with the same cycle structure {l"^^ ,2"^ ,3"^ , . . . , N"'' ) , are identical (Ap- 
pendix C). 

We will use the theory of group characters to decompose f' and 7' into basis functions which trans- 
form under irreducible representations of Sn- The theory of group characters is briefly discussed in Ap- 
pendix C. To determine the characters, Xj:^ ^'^ ^'■■■.Jv ") Qf ^j^g representation under which f' transforms we 
need to calculate how many of the elements f' transform into themselves under a permutation in the class 
(j^i^i 2''2,3''3, . . . , AT''"). As any number occurring in a cycle of length greater than one is transformed into 
another number, only those elements of f', rj, for which the index i occurs in a cycle of length one transform 
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into themselves. As there are vi cycles of length one, there will also be vi elements of r' which transform 
into themselves. Hence the character of the representation under which f' transforms under a permutation 
in the class (l^'S 2''^ 3''^ . . . , TV'' ) is given by 

Xr' -Vl- (69) 

The derivation of the characters, x^, ^'^ ^' • ■^ under which the 7' transforms is similar. For a given 
permutation in the class (I'^S 2"^^ , 3'^^ , . . . , N''^ ) all pairs of indices occurring in cycle of length one correspond 
to an element of 7', for which the indices i or j transform into themselves. There will be Vi{ui — l)/2 
such pairs. Also every pair of indices occurring in a cycle of length two correspond to an element of 7', 7^^-, 
for which the indices i and j transform into themselves. Thus the element transforms into itself. There are 
V2 such pairs. Thus the x^' ^'^ ^' ■ '^ ^\ under which the 7' transforms is the sum of these two numbers. 



'-■7' 

i.e. 

6.1.2 The Characters of the Irreducible Representations of Sn- 

Now let's turn to the characters of the irreducible representations. It can be shown[15] that for an element 
of an irreducible matrix representation which corresponds to a particular permutation with cycle structure 
(1^1 ,2''\'i"^ , . . . , N"^), the characters of the [N],[N -1, 1] and [N - 2, 2] representations are: 

(l-i,2'^2,...,Ar'^A') _ , 

X[iv-1, 1] > = Vi-l (71) 

(1-1,2-2,...,]^") _ (i^i - l)(z^i - 2) 

2] - 2 1- 2^2 - i • 

Armed with these characters, we can now see which irreducible representations of Sm f' and 7' span. 
In Appendix C we note that the characters of the representations under which f ' and 7' transform can be 
written as 

x{R)=Y.a^Xa{R), (72) 

a 

where R denotes the element of the group, x(^) is the character of the representation under which f' or 
7' transform, a labels the irreducible representation of 5jv , and Xa(-K) is the character of the irreducible 
representation denoted by a. The decomposition of Eq. (72) is unique where the coefBcient Ua is the number 
of times the irreducible representation labelled by a appears in the decomposition of the representation 
under which f' or 7' transform. 

Thus the f' vector must span a [N] and a [N — 1, 1] irreducible representation as from Eqs. (69) and 
(71) 

(l''i,2''2,...,JV''iV) _ (l''i,2''2,...,JV''iV) (l-i,2'^2,...,JV''A^) 

Xr' — ^ ^ X[N] T 1 X X[jv-1, 1] ' \''^) 
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i.e. 



a[N] = 1 (74) 
aiN-i, 1] = 1- (75) 

The 7' vector spans a [A''], a [A'' — 1, 1] and a [A'' — 2, 2] irreducible representation as from Eqs. (70) and 
(71) 

(l''i,2''2,...,JV''iV) _ (l-i,2-2,...,JV''iv) (1-1X',...,N-N) (l-i,2-2,...,JV''iv) 

A=y' — 1 X X[jv] i- 1 X X[jv-1, 1] i- 1 A X[Ar_2, 2] ' 

i.e. 

a[iv] = 1 (77) 

a[N-i, 1] = 1 (78) 

a[N-2, 2] = 1 • (79) 

Regarding the dimensionalities of the [N], [A'' — 1, 1] and [A'' — 2, 2] irreducible representations there are 
many rules to determine these. The best known is that the dimensionality of an irreducible representation 
of Sn is equal to the number of standard Young tableaux that the Young diagram for the irreducible 
representation allows. Perhaps the simplest rule however[17], is to fill the boxes of the Young diagram with 
numbers determined as follows. To the number of boxes in the row to the right of a box add the number 
of boxes in the column below the same box and then add one for the box itself. Once all of the boxes of 
the Young tableau have been filled with a number in this fashion, multiply all of the numbers in the boxes 
together. If we denote this product by Nt, then the dimensionality, d, of the representation is 

m 

d=— (80) 



from which we deduce that 



Nt 



, _ N\ _ N{N-3) 

^1 ~ [2m/{N{N-3))] ~ 2 ■ 



We note that 



d[,Y] +d[Ar_i, 1] = Af, (82) 
giving correctly the dimension of the f ' vector and that 

N{N-1) 

"[JV] + «[iV-l, 1] + "[JV-2, 2] = 2 ' ' ^ 

giving correctly the dimension of the 7' vector. 
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6.1.3 The Reduction of Df (R) by . 

Let Df' (i?) be the matrix representation of Sjq under which the f' vector transforms. The considerations of 
Sees. 6.1.1 and 6.1.2 then imply that Dfi{R) may be block diagonalized under a similarity transformation 
to 














(84) 



where D^^_]{R) is the one-dimensional irreducible matrix representation [A^] of Sn and Dg^~^' ^\R) is an 
{N — 1) X {N — l)-dimensional irreducible matrix representation [N — 1, 1] of Sn- 

If Wr' is a matrix which effects the reduction of Df>{R) to Eq. (84), then Wr' satisfies the equation 



Dr'{R) Wr^ [Wr'r'] = Dg,, {R) [W^'r'] . 



(85) 



Thus we have that the symmetry coordinates, Sr' , of the r' sector which transform under irreducible repre- 
sentations of the group Sn are given by 

Sf, = Wr'f' (86) 



and so writing 



(87) 



where W^^^ is a 1 x N dimensional matrix and W^^ ^' is axi {N — 1) x N dimensional matrix, then 



4^1 

^[JV-l, 1] 



(88) 



The symmetry coordinate column vector S^^^ is a one-element vector and transforms under £>^J (i?), the one- 
dimensional irreducible matrix representation \N\ of Sn- The symmetry coordinate column vector S^^i~^' 
is a (A/'— l)-element vector and transforms under £>^~^' (iJ), the (A^ — 1) x (A?^— l)-dimensional irreducible 
matrix representation \N —\, 1] of Sn- 



6.1.4 Reduction of D:^'(R) by W::^' 

. Likewise, if W~i is a matrix which effects the reduction of the representation, Diy'{R), under which the 7' 



vector transforms then 



where 



1%. D^, {R) W-,' [W^^-f'] = Ds^, {R) [W^'l'] 





( 





\ 


Ds^>{R) = 












I 








(89) 



(90) 
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£)^|(J?) is the one-dimensional irreducible matrix representation [A''] of ^jv, ^\R) is the (A'' — 

1) X (A'' — l)-dimensional irreducible matrix representation [N — 1, 1] of Sn , and ^k-^) the 

{N{N - 3)/2} X {N{N - 3)/2}-dimensional irreducible matrix representation [N - 2, 2] of Sn- 

Note that while D^,{R) and D^,{R) both belong to the a irreducible representation of Sn, where a is 
the partition [N] or [A'' — 1, 1], they are not necessarily the same matrices. Generally they are related by 
an equivalence transformation 

D^,{R)=tD^,{R)t-\ (91) 

where t is the transformation matrix. 

Thus we have that the symmetry coordinates, S^', of the 7' sector which transform under irreducible 
representations of the group Sn are given by 



S^' = W^, 7' 



(92) 



and so writing 



^if-i, 1] 

^_[W-2, 2] 
7 



(93) 



where is a 1 x A^(A^ — l)/2 dimensional matrix, is an (A^ — 1) x N(N — l)/2 dimensional 

matrix and W^, ^' is an N{N — 3)/2 x A^(A^ — l)/2 dimensional matrix, then 



■7 

gW-i, 1] 

glN-2, 2] 



(94) 



where 



N 



j=l i<j 



(95) 



The symmetry coordinate column vector Sl^^ is a one-element vector and transforms under £>^|(i?), the 
one-dimensional irreducible matrix representation [A''] of Sn, while the symmetry coordinate column vector 
S^~^' is an (A''— l)-element vector and transforms under -Dg^"^' {R), the (A''— 1) x (A''— l)-dimensional 
irreducible matrix representation [A'' — 1, 1] of Sn- The symmetry coordinate column vector S^~'^' 
an N{N - 3)/2-element vector and transforms under D^^_~'^' ^>(i?), the {N{N - 3)/2} x {N{N - 3)/2} 
dimensional irreducible matrix representation [N — 2, 2] of Sn- 



is 
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6.1.5 The Full Symmetry Coordinate Vector S. 

Symmetry coordinates S^, and S^, are said to belong to the same species when a = (3. Now in what follows 
we will find it useful to form a full symmetry coordinate vector as follows: 



P 




( \ 

g[N-2, 2] 



(96) 



where the orthogonal matrix 



/ 1 







10 



and 



^[JV-1, 1] 



10 
1 




1] 



(97) 



[N-1.. 1] 



and 



^[JV-2, 2] ^ 5i^-2. 2] 



(98) 



In S, symmetry coordinates of the same species are grouped together. 
6.2 Normal Coordinates 

So how do symmetry coordinates simplify the solution of Eqs. (46), (47), (48), (49) and (52) for the frequencies 
and normal modes? Symmetry coordinates have maximum utility when 

1. all transformation matrices t, which relate the equivalent irreducible representations of the symmetry 
group under which symmetry coordinates of the same species transform (see Eq. (91)), are equal to 
the identity matrix. The irreducible representations are not just equivalent, they are identical. 

2. the symmetry coordinates transform under unitary representations of the symmetry group (here Sn)- 
In our case we wish to have real symmetry coordinates and normal mode coordinates, and so we require 
that the symmetry coordinates transform under orthogonal representations of Sn- 

Now consider applying a transformation W to Eqs. (46), (47), (48) and (49), where 



W = P 














V' 





W-. 



N] 



2] 



(99) 
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where 

S = Wy'. (100) 

6.2.1 The Essential Equation and the Orthogonality of W. 

In Appendix E we show that condition 2 above implies that W has to be an orthogonal matrix, i.e. W 
satisfies 

WW'^ = W'^W = 1 . (101) 
In particular we show that W is an orthogonal matrix when 

W^,[WS,,V = 1^, (102) 

where X.' is r' or 7' and is the unit matrix. 

Thus Eq. (102) is the essential equation for W^i to satisfy. 

From Eq. (99) 

which with Eq. (101) gives us 

E [^x'fW^' = Ix' , (104) 

a 

where Ix' is the N x A'^-dimensional unit matrix when X' = f' or the N{N— 1) /2 x N{N— l)/2-dimensional 
unit matrix when X' = 7' . Prom Eq. (102) and the automatically satisfied Eq. (D4) we find 

[W^,fW^, m,fW^> = ^a,a' m^fW^, , (105) 

and so Eq. (104) is a decomposition of the identity into primitive idempotent projection operators which 
annul each other on both sides. The [W^^/]-^VF^, project out orthogonal irreducible subspaces while the 
are the Clebsch-Gordon coefficients of this decomposition. 

6.2.2 The Motion Associated with the Symmetry Coordinates about the Lewis Structure 
Configuration. 

Equation (104) can be written as 



Ix'=Y,Yl [(^^' )« ' ^^^"^^ ^ = 1 ^^^^ " = [^] - (106) 

1<^<N-1 when a = [N - 1, 1] , 

1 < ^ < N{N - 3)/2 when a=[N-2, 2], 



" ^ l<^<N-lwhena = [N-l, 1], 



21 



and the symbol ^ denotes the ^^"^ row of W^, . Thus we can write 

X' = Ix,X' = Y,[W^,V{W^,.X')^Y.Y. ((^x')«-^') [{W^'W, (107) 

where X' is f' or 7' of Eqs. (39) or (40) respectively. According to Eqs. (D4) and (102) the [(W^,)^]'^ are 
column vectors which form a complete set of orthonormal basis functions where, from Eqs. (D4) (102), one 
has 

{W^')d[Wi'Y\' = ■ (108) 

Prom Eqs. (88) and (94) 

[S'k']i = {W^')i-X' (109) 
and so Eq. (107) is a decomposition of the f' and 7' vectors in terms of the symmetry coordinates, i.e. 



la 



(110) 



where 



xT = [si']d{w^')ir (111) 

expresses the motion associated with the symmetry coordinate [S^,]^ in the original internal displacement 
coordinate, X' (f' or 7'). Equation (110) can be written concisely as 

y' = W^S, (112) 

where y' is defined in Eqs. (38), (39) and (40), S is given by Eqs. (96), (97) and (98) while W is given by 
Eq. (99). Equation (112) may be directly derived from Eq. (100) using the orthogonality of the W matrix 
(Eq. (101)). 

Using Eqs. (110) and (111) in Eqs. (22), (32) and (33) we find that 

/ 



y = 



where according to Eq. (Ill) 



r 
7 



(113) 



while 



and 



1 V l<i< N 



m']ij = i V i<i,j<N. 



(114) 
(115) 

(116) 

(117) 



Equations (113), (114), (115), (116) and (117) express the motion associated with the symmetry coordinate 
about the Lewis structure configuration. 
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6.3 The Reduction of the Eigensystem Equations in the Symmetry Coordinate 
Basis . 

6.3.1 The Central Theorem 

Let us start out by defining Qw to be the Q matrix in the symmetry coordinate basis, i.e. 

Qw = WQW^, (118) 

where Q is F, G or FG, and Qw is Fw, Gw or {FG)w ■ It has been shown elsewhere[18] that since F, 
G and FG are invariant matrices (under Sn), when items 1 and 2 above are satisfied then 



Qw = 



^ cr^N] ^ i[N] ^ 

1]«'I[N-1, 1] 

V '^?N-2, 2] ® I[N-2, 2] J 



(119) 



where the symbol (g) indicates the direct product, crj^j and are 2x2 dimensional matrices and 

""(AT— 2 2] ^ single element (a number). The matrix Ijn] is the [A^]-sector identity matrix, simply the 
number 1, I[n-i, i] is the [N — 1, l]-sector identity matrix, the {N — 1) x {N — l)-dimensional unit matrix 
and I[N-2, 2] is the [N — 2, 2]-sector identity matrix, the N{N — 3)/2 x A/'(A''— 3)/2-dimensional unit matrix. 
The matrix elements 

kg]x(,x^ = {W^Oi ^^i^i K^xO^f ' (120) 
where X[ and X2 are f' or 7' (only 7' for the [N — 2, 2] sector). The Q-matrix quadrant, Qx{x^, is given 
by Eqs. (59), (60), (61), (62) and (63). Note that although we are not summing over the repeated index 
^ in Eq. (120), the [cr^]x;,x^ are independent of the W^, row label, ^. If, when we calculate W^i , the 
matrix element [(t^]x[,x!^ turns out to depend on the W^, row label, ^ , then we know that we have made 
a mistake calculating W^, . This is a strong check on the correctness of our calculations in Sect. 7. 

6.3.2 The Reduction of the Eigenvalue Equation in the Symmetry Coordinate Beisis . 

Transforming the basic eigenvalue equation, Eq. (47), to the symmetry coordinate basis we have 

WFW'^ WGW'^ Wb = Fw Gw c^''^ = h c^^^ , (121) 

where 

Fw = WFW'^ Gw = WGW'^ c^'''> = Wb. (122) 
Under this transformation we can also write Eq. (46) for the b^^ normal-mode coordinate as 



[q% = b^W^Wy' = c^''^ -S, (123) 



T 



where [q']b is now directly expressed in terms of the symmetry coordinates. 
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Prom Eqs. (119) and (121) the normal-coordinate coefficient vector, c^''\ has the form 

— Oct, 1] ' J 5$) 



C 



(124) 



where the c" satisfy the eigenvalue equations 

cT^^c" = A„c- . (125) 

Note that cr^^ and cr^j^^i arc 2 x 2-dimcnsional matrices, while (t^^_2 2] ^'^ ^ onc-dimcnsional ma- 
trix. Thus there are five solutions to Eq. (125) which wc denote as {A^j, c[^'} , {A|^_j^ ""^^ and 
{A[7v-2. 2], c'-'^"^' ^1} , where X' (r' or 7' , only 7' for the [TV — 2, 2] sector) labels the rows of the elements 
of the column vector c" . The normal-coordinate label, 5, has been replaced by the labels a, ^ and ± on the 
right hand side of Eq. (124), while the elements of the column vectors 1^ are 

[n]v = kv^ (126) 

with l<(,,ri<N-l when a = [iV - 1, 1], or 1 < ^, r? < N{N - 3)/2 when a = [N - 2, 2] . The 
for the [A''] and [N — 1, 1] sectors determine the amount of angular-radial mixing between the symmetry 
coordinates in a normal coordinate of a particular a since from Eqs. (96), (98) and (123) 

W]b = [clh'[S$']i + [clW[S^.]i. (127) 

For the [A'' — 2, 2] sector we have 

= [cl^-^' ^1]^, [Sf-^- \ , (128) 

i.e. the [N — 2, 2] sector symmetry coordinates are the [N — 2, 2] sector normal coordinates up to a nor- 
malization constant, [cl^~^' ^Ij^y/ (see Sec. 6.3.3). 
If we write c" as 

c" = c" X c« (129) 

where c" is a vector satisfying the normalization condition 

[c5]'^c^=l (130) 

and c" is a normalization factor ensuring that Eq. (49) is satisfied (see Sec. 6.3.3 below), then the reduced 
eigenvalue equation, Eq. (125), determines c" alone. 
With the [N - 2, 2] sector, Eq. (130) yields 

c[J^^2l = 1 (131) 

and so Eq. (125) yields directly 

^^[^-2, 2] = 2] • (132) 
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As for the [N] and [A'' — 1, 1] sectors, if we write 



f a 

C± — 



COS 0% 

sin 612 



then from Eqs. (127), (129) and (133) 

[q']b = cl {cosei [S^,]^ + smBl [S^,]^) . 

Writing 



FG 



7 



then Eq. (125) may be written as 



008 6*2 



from which we derive 
A™ = 



Equations (136) and (137) then give us 



tan6iS 



6.3.3 The Normalization Condition. 

Prom Eqs. (49), (122), (119), (120) and (124) the c" also satisfy the normalization condition 



fc°l^a^c« = 1 . 



Thus Eqs. (129) and (139) mean that 



c = 



For the \N - 2, 2] sector, Eqs. (131) and (140) mean that 



2] 



'[iV-2, 2] 

As for the normalization constant, c%, of the [A''] and [A" — 1, 1] sectors, Eqs. (133) and (140) 

1 





T 




/ cos6»^ ^ 




{ COS 6*^ \ 


I sin - 




I sin ) 
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6.3.4 The Normal Coordinates. 



Thus the total transformation matrix, V , from the internal displacement coordinates to the normal coordi- 
nates is 

V = CW, (143) 

where W is given by Eq. (99) and 



C = 



' cos^f] 


cfsin^fl 




































^Icos^f-^' ^1 


4^-^' ^in^f-^' ^1 


































2] 



(144) 



and so the normal-coordinate vector, q' , is given by 



Q = 



V 9 



j'lN-1, 1] 
/[^-2, 2] 





[AT] 






[N-l, 


1] 




+ 








1] 


9 







,[JV-2, 2] 




cL^lcos^?fl5r 



+ 

[W] „[N] 



-[^Isin^^flsi^l 



cosr'^s. 



c_ sm I 



+ 
[N] 





cos 


9f-^' 


l]^[Af-l, 1] 


[N-l, 1] • 

+ sm 




1] ^^l^f-l. 


1] 




^' cos 






1 [N-l, 1] • 

+ sm 




1] 


1] 



4iV-2, 2]gr]A'-2, 2] 



(145) 



6.3.5 The Motion Associated with the Normal Coordinates about the Lewis Structure Con- 
figuration. 

Consider now the inverse transformation to express the internal displacement coordinates in terms of the 
normal coordinates. From Eq. (101) 



I = w'^w = w'^c-^cw 



(146) 



so that 



y' = W'^C-^ {Vy') = W^C-^ q' , 



(147) 
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where 
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and 



s{e°') = sin{e%-ei). 



(149) 



Thus from Eqs. (38), (147), (148) and (149) 
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Thus from Eqs. (22), (29), (32), (33) and (150) we find that 



(151) 



where 



(152) 
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with If I and 1-^' given in Eqs. (116) and (117) respectively. The vectors +yj", -J/j" and y^'^ ^' are 



s(6'«) 



for a = [N] or [iV- 1, 1] , and 



cos 6*° 
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s{0") 



■ cos 6[ 



(153) 
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(155) 



Equations (151), (152), (153), (154) and (155) express, in terms of the internal coordinates r and 7, the 
motion associated with the normal coordinates, q' , about the Lewis structure configuration . 



7 Determining the Symmetry Coordinates 
7.1 Primitive Irreducible Coordinates 

We determine the symmetry coordinates, and hence the WV and Wj of Eqs. (86) and (92) respectively, in a 
two-step process: 

a) . Determine two sets of linear combinations of the elements of coordinate vector f' which transform 

under particular non-orthogonal [N] and [N — 1, 1] irreducible representations of ^jv • Using 
these two sets of coordinates we then determine two sets of linear combinations of the elements 
of coordinate vector 7' which transform under exactly the same non-orthogonal [N] and [A''— 1, 1] 
irreducible representations of Sn as the coordinate sets in the f' sector. In this way we satisfy 
item 1 of Subsec. 6.2 above. Then another set of linear combinations of the elements of coordinate 
vector 7' which transforms under a particular non-orthogonal [A''— 2, 2] irreducible representation 
of Sn- These sets of coordinates, which we term primitive irreducible coordinates for reasons 
discussed below, have W matrices which automatically satisfy Eq. (D4). 

b) . Appropriate linear combinations within each coordinate set from item a), above are taken so 

that the results transform under orthogonal irreducible representations of Sn ■ These are then 
the symmetry coordinates of Section 6 above. Care is taken to ensure that the transformation 
from the coordinates which transform under the non-orthogonal irreducible representations of 
Sn of item a), above, to the symmetry coordinates which transform from the under orthogonal 
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irreducible representations of Sn preserve the identity of equivalent representations in the f' and 
7' sectors to ensure that item 1 of Subsec. 6.2 above continues to be satisfied. In this way both 
items 1 and 2 of Subsec. 6.2 above are satisfied. 

What are the advantages of this two-step process? In step a), above we set out to find sets of coordinates 
transforming irreducibly under Sn which have the simplest functional form possible. This is why we call them 
primitive irreducible coordinates. In step b). above, we find appropriate linear combinations of the primitive 
irreducible coordinates so that they transform under orthogonal representations of Sn- Thus the symmetry 
coordinates are composed of the building blocks of the primitive irreducible coordinates. Furthermore, we 
choose one of the symmetry coordinates to be just a single primitive irreducible coordinate, and so it describes 
the simplest motion possible under the requirement that it transforms irreducibly under Sn- The succeeding 
symmetry coordinate is then chosen to be composed of two primitive invariant coordinates and so on. In 
this way the complexity of the motions described by the symmetry coordinates is kept to a minimum and 
only builds up slowly as more symmetry coordinates of a given species are considered. 

7.2 The Primitive Irreducible Coordinate Vector, Sr • 

The primitive irreducible coordinates, Sri, of the f' sector which transform under irreducible, though non- 
orthogonal, representations of the group ^jv are given by 

Sri=Wrif', (156) 

where 

= J-^_,, , (157) 

\ ^r' J 

W^^^ is a 1 X A'' dimensional matrix and W^^ ^' is an {N — 1) x N dimensional matrix. Hence we can 
write 

Since Wr> is a matrix which cfi^ccts the reduction of the representation, Dr'{R), under which f' trans- 
forms, to particular irreducible, though non-orthogonal, representations of the group ^jv under which Sr' 
transforms, then 

Wr' Df' {R) Wf! Sf' = Dg ^ (J?) 'Srs , (159) 



where 











(160) 



In the above, -D^J {R) is the one-dimensional irreducible matrix representation [N\ of Sn and ^' ^\r) 
is an {N — 1) X {N — l)-dimensional, non-orthogonal, irreducible matrix representation [N — 1, 1] of Sn- 
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Thus the primitive irreducible coordinate column vector Sf^, is a one-element vector and transforms under 
£>^J(i?), a one- dimensional non-orthogonal irreducible matrix representation [N] of S^. The symmetry 
coordinate column vector ^' is a (A'' — l)-element vector and transforms under D^_^ ^' ^\r), the 
aforementioned {N — l)x {N— l)-dimensional, non-orthogonal, irreducible matrix representation [N — 1, 1] 
of ^JV. 



7.3 The Primitive Irreducible Coordinate Vector, S- 



Likewise the primitive irreducible coordinates, S~', of the 7' sector which transform under irreducible, 



though non-orthogonal, representations of the group Sn are given by 

S^, = W^' 7' 

and so writing 



rv y 



(161) 



(162) 



where VK^' is a 1 X N{N — l)/2 dimensional matrix, ^' is an (A^ — 1) x N{N — l)/2 dimensional 

matrix and W^^ ' is an N{N — 3)/2 x N{N — l)/2 dimensional matrix, then 



Sy = 



W^-i, 1] 



7 



where 



N 



(163) 



(164) 



Since Wy' is a matrix which effects the reduction of the representation, D^'{R), under which the 7' 
vector transforms, to particular irreducible, though non-orthogonal, representations of the group ^jv under 
which Sy' transforms, then 

W^' (R) W^} S^> = Dg_^ {R) S^> , (165) 



where 










\ 







dI;^' \r) 






[ 





D^- \R) ) 



(166) 



In the above, D^_^ (R) is the one-dimensional irreducible matrix representation [N] of Spf, ^ ^\R) 

is an (N — 1) x (N — l)-dimensional, non-orthogonal, irreducible matrix representation [N — 1, 1] of Sn 

[A^ — 2, 2] 

and D-a_ ' (i?) is an {N{N — 3)/2} x {N{N — 3)/2}-dimensional, non-orthogonal, irreducible matrix 
representation [A^ — 2, 2] of Sn- 
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Note that according to item a), above, 



(167) 



where a is the partition [N] or [A'' — 1, 1] , 



Thus primitive irreducible coordinate column vector S-^i is a one-element vector and transforms under 

fjV] 

D (R), the one-dimensional irreducible matrix representation [A''] of Sn, while the primitive irreducible 
coordinate column vector ^' is an (A'' — l)-element vector and transforms under D^^ ^' ^\r), an 
{N — 1) X {N — l)-dimensional, non-orthogonal, irreducible matrix representation [N — 1, 1] of Sn- The 

[JV— 2, 2] 

primitive irreducible coordinate column vector S^, ' is an N{N — 3)/2-element vector and transforms 

[AT— 2 21 

under D' ' '{R), the {N{N - 3)/2} x {N{N - 3)/2}-dimensional, non-orthogonal, irreducible matrix 
representation [N — 2, 2] of Sn- 

We also note that according to Eq. (D4) 



7.4 The Full Primitive Irreducible Coordinate Vector, S 



(168) 



Like for the discussion of the symmetry coordinates in Subsection 6.1 above, it is useful to form a full 
primitive irreducible coordinate vector that groups primitive irreducible coordinates of the same species 
together as follows: 



S = P 




iJfif 



-^[N-l, 1] 
1] 

-^[^-2, 2] 



-^[N-l, 1] 



v 



2] 



(169) 



where the orthogonal matrix P is given by Eq. (97) and 



-AN] 



1] 



and 



(170) 



Now consider applying a transformation W to Eqs. (46), (47), (48) and (49), where 



W 



Wf' 












(171) 



where 



Wy'. 



(172) 
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Eq. (172) implies that the representation under which the displacement coordinate vector S transforms is 



Ds = WDy,{R)W 



(173) 



where Dyi{R) is given by Eq. (El). The non-orthogonal matrix representation D-g[R) is a block diagonal 
matrix of irreducible, though non-orthogonal, representations of the form 
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V 

where we have also used Eq. (167) from above in the last step. 



(174) 



7.5 The Transformation, U , from Primitive Irreducible Coordinates to Sym- 
metry Coordinates. 

The symmetry coordinate vector S will be related to the primitive irreducible coordinate vector S' by a 
non-orthogonal linear transformation U, i.e. 

S = US (175) 

and likewise Ds{R) will be related to D-g{R) by 

Ds{R) = UDgiR) U-^ . (176) 



Thus from Eqs. (100), (172), (E2) and (173), W will be related to W by the same transformation, i.e. 



W = UW , 



(177) 



where 
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Tj[N-2, 2] 



(178) 



so that 



W = 



..[jv-i, 1] yr?[N-h 1] 



jj\N-1, 2] TT7[f^-2, 2] 



(179) 
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Note that from Eqs. (E8), (174) and (176), 

U2, = U^', , (180) 

where a is [N] or [iV — 1, 1] and Afj is a number (see below). 
Comparing Eq. (179) with Eq. (99) we see that 

W^' = U^' (181) 

and according to Eq. (102) J7^, must satisfy 

t/x' {wUWx'f} [U^.f = Ia- (182) 

Thus to determine C/^, we need to know beforehand. Note that when we use Eq. (180) in 

Eq. (182) we must have 



when a is [N] or [AT - 1, 1] 



7.6 Symmetry Coordinates Belonging to the [N] Species 
7.6.1 The f' Sector 

The Primitive Irreducible Coordinate The primitive irreducible coordinate Sfi which is invariant 
under Sn is readily seen to be 

N 



s':!'=Y.fl (184) 

i=l 



[AT] 

From Eq. (158) we can identify W^,, as 



[W[V], = [U,]i (185) 

where lf> is defined in Eq. (116). We leave the calculation of s]^ ^' ^' and W^^ ^' ^' to a subsequent 
paper [10]. 

The Symmetry Coordinate One readily sees that 



W^^\W^^Y = N (186) 



and so from Eq. (182) 

ttW - 



[7^1 = 4?- (187) 



Thus from Eqs. (175), (181), (184), (185) and (187) 



[W^fl], = -^[V], (188) 



and 

N 



We leave the calculation of wl.^ ^' and ^' to a subsequent paper [10]. 
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Motions Associated with Symmetry Coordinate . According to Eqs. (113), (116), (114) and 
(188) the motions associated with symmetry coordinate S^^^ in the unsealed internal displacement coordi- 
nates r about the unsealed Lewis structure configuration 

Too = D^Uhof^lr' (190) 

are given by 

rl^l = D'/' f'l^l = D^/' S^^^ [(W^jf 1)]^ = a,, ^^4^1 1,, . (191) 

Thus wc sec that the motions associated with symmetry coordinate S^!^^ involve symmetric-stretch motions 
about the Lewis structure configuration, where all of the radii expand and contract together. 
We leave the calculation of r^^~^' to a subsequent paper [10]. 

7.6.2 The 7' Sector 

The Primitive Irreducible Coordinate According to item a), above, S^y/ should transform under 
exactly the same non-orthogonal irreducible [N] representation of Sn as s]^^ . Since S^^^ is invariant under 
all of the elements, we simply require that S:=^, be invariant under ■ 
Now according to Eq. (4) 

7ij =ri.fj, (192) 

where fi is the unit vector from the center of the confining field to particle i. Now J2f=i invariant 
under Sn and so {J2iLi ^i) • (^jLi r?) invariant under Sn ■ The contribution from the terms with i = j, 
i.e. ri.ri = N , is also invariant on its own. Since this is just equal to the number N, we subtract 

this term from {J2iLi ^j) • (Sj^i '^O ^^'^ consider the invariant quantity {J2iLi '^i) • i^jj^i r?) • According 
to Eq. (192) this equals J2iLi ^j^i lij = 2 Yl!j=2 ^i<j 7^ ^^'^ npon using Eq. (33) we obtain [N{N - 
7oo +^^^^ Sj=2 Si<j ^'ij ■ Siiice the first term is invariant under Sn, J2jL2 Si<j ^'ij invariant under 
Sn, and so we identify the primitive irreducible coordinate of the 7' sector which transforms under the [N] 
representation as 

j=2 i<j 

From Eqs. (40), (163), (164) and (193) 

[Wl^l^hj = , (194) 

where 1:^' is defined in Eq. (117). We leave the calculation of ^\ Wl^~^' andW^^"^' 

to a subsequent paper [10]. 

The Symmetry Coordinate One readily sees that 

Wl^^ [W^']- = f ^ [Wl^%, [W^'], = , (195) 
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where we have used Eqs. (194) and (117) in the last step, and so from Eq. (182) 



Thus from Eqs. (175), (181), (193), (194) and (196) 

AN] 



and 

AN] I 2 



We leave the calculation of ^\ ^' sl^~^' and S'i^"^' ^'to a subsequent paper [10]. 



Motions Associated with Symmetry Coordinate . According to Eqs. (113), (115), (117) and 
(197) the motions associated with symmetry coordinate in the unsealed internal displacement coordi- 
nates 7 about the Lewis structure configuration 

7oo=7oolr (199) 

are given by 

^[N] _ J_ V[iv] _ J_ a^N] UurlNU^T _ / ^ c[^] i_, (200) 

^ ~ VD^ ~ VD ^' ^' ~\l N{N -1)D • ^^""^ 

Thus we see that the motions associated with symmetry coordinate Sl^^ involve symmetric-bend motions 
about the Lewis structure configuration, where all of the angles open out and contract together like a parasol, 
but with the radii unchanged. 

We leave the calculation of 71^^!' ^1 and 'y^^~'^' ^1 to a subsequent paper[10]. 

8 The Frequencies and Normal-Mode Coordinates of the 
[N] Species of the System. 

8.1 The G and FG matrices in the Symmetry-Coordinate Basis. 

We can use the W^} matrices of Eqs. (188) and (197) to calculate the reduced G and FG matrix elements, 
[f ^]]x(,x^ and [<t^]xj,x^ , of Eq. (120). An outline of this calculation is to be found in Appendix F, 
with the following results. 

8.1.1 The Matrix Elements [CT^]]xi,x^ • 

Using Eqs. (59), (60), (61), (62), (63), (66), (67), (68), (120), (188) and (197) we derive 



(201) 
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We leave the calculation of the matrix elements and <t[j^_2, 2] ~ [""^^2, 2]W,y' *° ^ 

subsequent paper [10]. 

8.1.2 The Matrix Elements [fT^]xi,xi • 

Using Eqs. (59), (60), (61), (62), (63), (64), (65), (120), (135), (188) and (197) we derive 



^ ['^fi^]h',r' = {a + bN) [^fr^]]r',r = V2(7V-1) + f /) ^ 

^ Hi^]W,r' = VW^) + jd^ ['r[r^]W,^' =[9 + 2{N - l)h + i^i^I^ ?^ 



(202) 

We leave the calculation of the matrix elements i]]x(,x^ and cr^_2 2] = 2]]t.T ^'^ ^ 

subsequent paper [10]. 

8.2 The Frequencies and Normal Modes. 

Using Eqs. (202) in Eq. (137) we obtain Aj^j . The frequencies are then determined from Eq. (52). The f'-^' 
mixing angles. (9[^1 for the [N] species arc determined from Eq. (138). The normalization constant c'^l of the 
reduced normal-coordinate coefficient vector, c'^l, of Eqs. (124) and (129) are determined from Eqs. (141) 
and (142). One then determines the normal mode vector, q' , through Eqs. (145), (189) and (198). Thus we 
arrive at 

, 1 f;^. + i"Ui„,(M f^j2% (m 

^ i=l V ^ ' j=2 i<j 

We leave the calculation of q'^^ ^' ^\ q'^^~^' and q''^^^' ^' to a subsequent paper[10]. 

8.3 The Motions Associated with the Normal Modes. 

From Eqs. (151), (152), (153), (154), (155), (188) and (197) 



where yoo is given by Eq. (152). The ^ sum for the [N] species only includes one term and so 



ID^ sin 6* 
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Clho 



D3 sine^^^ 




[JV] 



y_ J-r-' ) 



-7 



[JV] 



cos^l 



g_ 1- 



]]N{N-1)D s(0[iv])cL^l 
We leave the calculation of rT"[^~^' , Tjf^~^' and 7j^~^' ^' to a subsequent paper[10]. 



(208) 
(209) 



9 Summary and Conclusions 

In earlier papers[2, 4, 8] we have developed the method of dimensional perturbation theory[3] at low or- 
ders to examine the energies of quantum-confined systems such as atoms, quantum dots and Bose-Einstein 
condensates, both the ground and excited states. Dimensional perturbation theory (DPT) has many advan- 
tages. These include the fact that the number of particles, N , enters into the theory as a parameter, and 
so it is easy to calculate results for an arbitrary number of particles. Also the theory is a beyond-mean-field 
method, directly accounting for each interaction, rather than some average representation of the interactions. 
This means that it is appropriate for confined systems of both weakly interacting and strongly interacting 
particles. In the case of a trapped gaseous atomic Bose-Einstein condensate (BEC) the system is typically 
a weakly interacting system for which a mean-field approximation is valid. However, if A'' is increased, or 
Feshbach resonances are exploited to increase the effective scattering length, or the system is squeezed in 
one or more directions, it will transition into a strongly-interacting regime where the mean-field approach 
breaks down. Such systems have been created in the laboratory. Dimensional perturbation theory, how- 
ever, is equally applicable to both regimes and may be used to study the transition from weak to strongly 
interacting regimes. Even at lowest order DPT includes beyond-mean-field effects. 

In this paper we have begun to address the problem of calculating the zeroth-order DPT wave function for 
quantum-confined systems, and it is a considerably expanded discussion of work briefly outlined in a previous 
letter [9]. At large dimensions the Jacobian- weighted wave function of the system becomes harmonic, with 
the system oscillating about a configuration termed the Lewis structure. Notwithstanding the relatively 
simple form of the large-dimension, zeroth-order wave function, it includes beyond-mean-field effects. 

With the wave function many more properties of the system become accessible. The normal mode 
coordinates of the zeroth-order wave function reveals the nature of the excitations of the system and the 
lowest-order wave function yields expectation values and transition matrix elements to zeroth order in DPT. 
For macroscopic quantum-confined systems, such as the BEC, the wave function is uncloaked in an explicit 
way as the density profile may be viewed in a direct fashion. This is readily calculable at zeroth order in 
DPT from the zeroth-order wave function. Also, calculating energies and wave functions to higher orders in 
1/D requires as input the lowest-order wave function. 

In this paper we limit ourselves to setting up the general theory of the zeroth order DPT wave function 
for a quantum-confined system. We illustrate this by applying it to the breathing and center of mass modes 
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of the [N] species. In the case of the BEC the ground state energies and excitation frequencies of these 
modes have been calculated using DPT in a previous paper. The modes of the [N — 1, 1] and [A'' — 2, 2] 
species are derived in a subsequent paper. The presence of these additional normal- mode coordinates in the 
theory indicates the existence of modes of excitation of the BEC beyond those which are typically considered. 
The zeroth-order density profile of the ground state condensate is to be considered in another paper, and is 
relatively simple to calculate once the zeroth-order wave function is known. 

This paper restricts itself to consideration of quantum-confined systems with a confining potential of 
spherical symmetry. The extension of this method to systems with cylindrical symmetry is relatively straight- 
forward and is to be discussed in subsequent papers. This is particularly important for the BEC as almost 
all current laboratory traps have cylindrical symmetry. 

As we have discussed in an earlier paper [4] considering energies and excitation frequencies of the BEC, 
low-order DPT calculations of confined interacting particles, and this includes strongly interacting systems, 
are accurate out to N equals a few thousand particles. To obtain more accurate energies, and to extend 
the calculation beyond a few thousand particles we need to go beyond these low-order calculations of the 
energies and excitation frequencies. These calculations are facilitated by the present paper since, as we have 
already mentioned, it requires the zeroth-order DPT wave function. This will be pursued in later papers. 
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A The Wilson FG matrix method 

In this appendix, wc derive the Wilson FG matrix method[14] which is at the heart of our obtaining the 
normal mode frequencies, the first-order energy correction, the normal coordinates and the wave function at 
lowest order. The derivation involves a transformation to the set of coordinates called normal-coordinates 
in which both the differential term and the potential term of Eq. (45) are diagonal. We begin by defining 
a symmetric transformation, A, that transforms from the vector y', defined by Eq. (38), to z' (z' = Ay'). 
A is an active transformation and it has the property that it diagonalizes the symmetric G matrix to unity. 
That is, A satisfies 



A^GA = I 



(Al) 



which can also be written as 



G = A-\A-Y, 



(A2) 
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where we have used the property that A is symmetric (A = A^) in the derivation of Eq. (A2). Hi then 
becomes 

Hi - —did,, + y^{A-YFA-^z'. (A3) 

Next we focus our attention on the potential term, the term involving the matrix F. We introduce 
another transformation, U {q' = Uz'), that diagonalizes the potential term while simultaneously leaving 
the differential term unchanged. This orthogonal transformation (U^U = I, where I is the identity matrix) 
leaves the differential term in the same form as in Eq. (A3), and the potential term becomes 

U(A-^)^FA-iU^ = A, (A4) 

where A is a diagonal matrix. That is. 

Hi - -\dl'd^' + ^g'^V- (A5) 

Prom Eq. (A5) it is clear that A is the diagonal matrix of the squares of the frequencies of the normal modes 
(see Eq. (52)). 

Now if we introduce the projection operators 

Pb = b"b"^, (A6) 

where the c"' element of the column vector b" is 

(&")c = She (A7) 

with b and c running from 1 to {N + l)/2, i.e. all of the elements of the column vector are zero except the 
6*'' element which is equal to unity, then 

A b"b"^ q' = [A],, b"[q% = X, b"[q']t, , (A8) 

where [q']i, , the b*"^ normal- mode coordinate, is given by 

W]b = b"^q' (A9) 

Since [q']b is just a multiplier on all sides of Eq. (A8), Eq. (A8) reduces to 

A6" = A6&", (AlO) 

Using Eq. (A4), Eq. (AlO) reads 

U(A-i)^FA-iU^ b" = Xb b" . (All) 
There are two important ways we can proceed from here. 
1 . If we now define 

b" = \5AB (A12) 
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then Eq. (A9), the equation for the 6*^ normal- mode coordinate, reads 

[q'], = B^AV = fi^G-V, (A13) 

where we have used Eqs. (A2) and (A4), and from Eq. (All) we find that B satisfies the eigenvalue 
equation 

GFB = \bB. (A14) 
The normalization of B is determined by the normalization condition for b", 

b"^b" = l (A15) 

(see Eq. (A7)). Using Eqs. (A2) and (A4) in Eq. (A15) we obtain 

B'^G-^B = 1. (A16) 

Notice how the constant matrix is the contravariant metric tensor for the space. This is to be 
expected from the form of Eq. (45). 

2. All of the work applying DPT to large--/V quantum-confined systems up till now has been focussed 
on calculating energies and adopts the procedure outlined in item 1. The frequency squared, Aft, is 
calculated from Eq. (A14) by solving 

det(GF- AB) =0, (A17) 

Since we are also interested in large- Z) wave functions and normal coordinates, our life would be easier 
if B were the coefiicient vector of the displacement coordinates contributing to the normal coordinate 
[q']b . Looking at Eq. (A13) we see that this is not so as there is a G~^ operation involved. Thus 
instead of following the procedure outlined in item 1 above, in this paper we solve directly for the 
coefficient vector, b, of the displacement coordinates contributing to the normal coordinate [q']b , 

b = G"^B. (A18) 
Thus in terms of 6, Eq. (A12) for the normal mode coordinate takes on the desired form 

W]b = b'y'. (A19) 
Using Eq. (A18) in Eq. (A14) we find that the coefiicient vector satisfies the eigenvalue equation 

FG6 = A6 6 (A20) 

with the resultant secular equation 

det(FG- AI) =0. (A21) 
Prom Eq. (A16) we see that the coefficient vector satisfies the normalization condition 

b^Gb=l. (A22) 
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B Gramian Determinants 



The Gramian determinant [19] is defined as: 



r = det(7ij), 



(Bl) 



where 7^ = • rj/riTj^ the angle cosines between the particle radii and r^, represents the elements of an 
N X N matrix. A related quantity used in the main text is the principle minor of the Gramian, r("\ defined 
as the determinant of the jij matrix with the a*^ row and column removed. 

A most challenging part of calculating the largc-D minimum and the F matrix elements in the systems 
discussed in this paper is handling the Gramian determinants and their derivatives. What makes these 
calculations feasible is the very high symmetry of the infinite-dimension, symmetric minimum. In this ap- 
pendix we discuss the Gramian determinant and show how to obtain its derivatives at the infinite-dimension, 
symmetric minimum. We will now demonstrate how this is done and summarize the results. 

In general = 1 and at the infinite-!? symmetric minimum all of the remaining direction cosines are 
equal, 7^- = 7;^. Hence, the Gramian determinant is an {N x N) matrix of the form 



= det 



7o 



/ 1 

7cx> 

V Too • • • ' 

in Appendices B and D of Ref. [2] we show that (Eq. (D4)) 



7o 



"Too 

1 / 



(B2) 



[1 + (JV-1)^J(1-^^) 



JV-l 



(B3) 

is simply related to the corre- 



The principle minor evaluated at the infinite-Z) symmetric minimum F^") 

00 

sponding Gramian determinant (B3) by A/" —> A/' — 1. 

To calculate the infinite-D symmetric minimum of the Gramian derivatives, we expand F in terms of its 
cofactors. The cofactor, denoted by Cij, of the element 7^ in F is (—1)*+-' multiplied by the determinant of 
the matrix obtained by deleting the i*^ row and j*'* column of F. We may then write F as 



N 

F = ^-f^jCij. 

Then the partial derivative of F in terms of the cofactor is 

-2C - 

Prom this equation, the partial derivative of F evaluated at the infinite- 2) symmetric minimum is 

dT 



= _2C7(f-i), 



(B4) 



(B5) 



(B6) 
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where we have defined the following determinant of an (A'' — 1) x (A'' — 1) matrix: 

\ 



= det 



/ 

/ 7 


Too 


Too 


"Too 


1 


Too 


"Too 




1 


V 7oo 







7o 



■• 7oc 

7<x> 1 



(B7) 



/ 



and the superscript {N — 1) simply indicates the size of the matrix. Prom this matrix one can show that the 
following recursion relation holds: 



or equivalently 



(Ar-2) 



- (iV - 2)C^^- 



(N-l) 



where the (A'') superscript in the notation F 



(JV) 



(B8) 
(B9) 

From the recursion 



again refers to the size of the matrix F 
relation (B9), one can easily prove by induction the conjecture that = ^^{1 — 7oo)^~^ so 

dT 



-27oo(l-7oo) 



N-2 



(BIO) 



The derivative of the principle minor evaluated at the infinite-D symmetric minimum is simply related to 
the corresponding Gramian determinant derivative (BIO) by N ^ N — 1. 

To summarize the above results, the following expressions are needed when calculating the minimum of 
the effective potential (26 and 27): 



or 

dlij 



-27oo(l-7oo) 



™ \N-2 



or'- 



-27oo(l-7oo) 



rU = [1 + (A^ - l)7oo](l - loor-' r(«)| =[1 + {N- 2)7^](l - 7oo) 



7^ \N-2 



(Bll) 



And when evaluating the F matrix elements at the infinite- 1? symmetric minimum, the following six 
second-order derivatives of the Gramian determinants are needed: 



= ^ 



= 



97« 



32p(o) 



97| 



-2(1 - 700^-^(1 + 3)7 
: -2(1 -700^-4(1 + (TV -4)7, 



= 27oo(l-7oo) 



7^ \N-3 



(B12) 



= 27cx>(l-7oo) 



TT, \N-A 



C The Symmetric Group and the Theory of Group Characters 

A group of transformations is a set of transformations which satisfy the composition law ah = c where a 
and h arc any two elements of the group and c also belongs to the group. A group also contains the identity 
element / such that al — la = a and for every element in the group its inverse is also to be found in the 
group. 
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The symmetric group Sn is the group of all permutations of N objects and, as such, has N\ elements. 
A permutation may be written as 

^12345678 
^23154768 

This denotes the following transformation 



(CI) 



Object before 
transformation 



is transformed 
to object 



1 
2 
3 
4 
5 
6 
7 



2 
3 
1 
5 
4 
7 
6 



A cycle is a particular kind of permutation where the object labels are permuted into each other cyclically. 
For example the cycle {abc) means that object a is transformed into object b, object b is transformed into 
object c and object c is transformed into object a. The cycle {abc) is termed a 3-cycle since it cycles three 
objects. Like wise (3479) is a 4-cycle and (5) is a 1-cycle, the latter transforms object five into itself. All A''! 
permutations of the group Sn may be decomposed into cycles. For example, the permutation of Eq. (CI) 
may be decomposed into cycles as 



12345678 
23154768 



= (123)(45)(67)(8) 



(C2) 



This consists of one 1-cycle, two 2-cycles and one 3-cycle. We can denote the cycle structure of a permutation 
by the symbol (f^^, 2"^^, 3"^, . . . , A''''"),[15] where the notation j"^^ means a cycle of length j and Uj equals 
the number of cycles of length j in that permutation. In the case of the permutation of Eq. (C2), its cycle 
structure is {l'^,2'^,3^). 

A matrix representation of a group is a set of nonsingular matrices, including the unit matrix, which 
has the same composition law as the elements of the group. The character of an element of an matrix 
representation of a group is the trace of the matrix. The character, as the trace of a matrix, is invariant 
under a similarity transformation and thus elements of equivalent representations have the same character. 
The set of all the distinct characters of the elements of an irreducible representation of the group uniquely 
specify the irreducible representation up to an equivalence transformation. The characters of irreducible 
representations of a group are termed simple characters. All elements of a group which are related by a 
similarity transformation are said to belong to the same class. The character of the elements of a group 
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belonging to a particular class all have the same character. Thus there are as many distinct characters for 
a group as there are classes. 

For the group Sn all elements with the same cycle structure belong to the same class and so all elements 
of a matrix representation of Sn with the same cycle structure have the same character. 

A reducible matrix representation of a group may be bought to block diagonal form by a similarity 
transformation, where the individual blocks are irreducible matrix representations of the same group with 
lower dimensionality. Thus the characters, x(^)) of a reducible group are the sums of the characters of the 
irreducible matrix representations into which it can be decomposed, i.e. 

x{R) = J2xp{R), (C3) 

p 

where R denotes the element of the group, p labels all of the irreducible blocks into which the reducible matrix 
representation of the group may be decomposed and Xp{R) is the character of the irreducible representation 
of the block labelled by p. Now in a particular reducible representation a given irreducible representation 
may be repeated along the diagonal aa times and so Eq. (C3) may be rewritten as 

x(i?) = ^a„Xa(J?). (C4) 

a 

The decomposition of x(^) of Eq. (C4) into simple characters Xa{R) is unique, i.e. there is not another 
decomposition of the form 

X(i?) = ^6oXa(i?), (C5) 
a 

where at least one of the ba is different from the corresponding a^. This follows from the fact that quite 
generally 

ac. = lY,xUR)x{R) (C6) 
and that the simple characters satisfy the orthogonality condition 

^xUR)X0{R) = h6af3, (C7) 

R 

where h is the number of elements in the group. For Eqs. (C5), (C6) and (C7) to be consistent we must have 

6„ = a„ . (C8) 

The irreducible matrix representations of may bo labelled by a Young diagram ( = Young pattern = 
Young shape). A Young diagram is a is a set of N adjacent squares such that the row below a given row of 
squares is equal to or shorter in length. The set of all Young diagrams that can be formed from N boxes of 
all possible irreducible matrix representations of Sn. 

A given Young diagram may be denoted by a partition. A partition, [Ai, A2, A3, . . . , Xn] is a series 
of TV numbers Ai such that Ai > A2 > A3 > . . . > Aat such that Ai + A2 + A3 + • • • + Aat = A^. The number 
Aj is the number of boxes in row i of the Young diagram. Thus the set of all possible partitions of length 
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N labels all of the possible irreducible matrix representations of Sn and so the irreducible representation 
labels a and f3 in Eqs. (C4), (C5), (C6), (C7) and (C8) above, for the group Sn may be taken to run over 
the set of all possible partitions. 

A Young diagram can have up to N rows. For an N-iow (one-column) Young diagram all of the AjS are 
non zero. However only one Young diagram will have N row; all of the rest will have less than N rows. Thus 
in all but one case the last few AjS will be zero. It is standard practice to drop the zeros and use a partition 
with less than A'' numbers. Thus the partition [3, 0, 0] labelling an irreducible representation of Ss is usually 
abbreviated to just [3]. 

D The Inequivalence of Subspaces a and cx' under Sn when ol^ ex! 
— A Useful Equation. 

Here we wish to note something that proves useful in Sec. 6.2 above. Consider the quantity 

W'^,T\wi,f Di^XR) , (Dl) 

where X' is f' or 7', a and o! denote the partitions [iV], [A'' — 1, 1] or [A^ — 2, 2] , and T is an unspecified 
matrix. The matrix W^,T[W^,]^ effects a similarity mapping from the irreducible a' and Sx' space onto 
an irreducible space denoted by a. Thus one can construct an irreducible matrix representation D^{R) 
belonging to the set of equivalent representations denoted by a such that 

D^{R) w^,r[w^,f = w^,r[w^,f d^^, {r) (D2) 

for all elements R of the group. According to Schur's Lemma however, any such similarity mapping, 
W^,T[VF^',]^, that satisfies Eq. (D2) is identically zero unless a = a', i.e. 

W^,r[W^,f = 0, V a^a'. (D3) 

Equation (D3) is a statement of the inequivalence of the irreducible spaces labelled by a and a' under Sn 
when a ^ a'. In particular we have 

W^,[W^,f = 0, V aj^a'. (D4) 
Schur's Lemma also implies -D" {R) ^ -Dg^, (-R) unless W^x'l^x'l"^ ^a', where /„/ is the unit matrix. 

E Proof of the Orthogonality of W . 

We shall show that condition 2 of See 6.2 above implies that W is an orthogonal matrix. This can be 
seen as follows. The displacement vector y' of Eq. (38) transforms under an orthogonal, though reducible. 
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representation of ^jv since y y' is invariant under Sn- Writing this orthogonal representation as 



Dy.{R) = 



Dr'iR) 











(El) 



Eq. (100) imphes that the representation under which the displacement coordinate vector S transforms is 

Ds = WDy,{R)W-^ . (E2) 
Since Dy' (R) is an orthogonal matrix representations of Sn then from Eq. (E2) we have 

D^,{R)Dy,{R) = W^D^{R)[W-YW-^Ds{R)W = 1. (E3) 
Prom Eq. (E3) we can write 

D'^{R)[WW^]-^Ds{R) = [WW^]-^ = D^{R)Ds{R)[WW^]-^ , (E4) 
where we have used the orthogonality of Ds{R) in the last step. Left multiplying by [Dg]~^{R) we obtain 

[WW'^] - ^ Ds {R) = Ds {R) [WW'^] " S (E5) 

from which we derive 

WW'^ Ds{R) = Ds{R)WW'^ . (E6) 

We notice that 

WW^ = I (E7) 

satisfies Eq. (E6). We will now show that it is the only solution. We note that Ds{R) is a block diagonal 
matrix of irreducible representations of the form 

Ds{R) = P 














£)[^l(i?) 

















D^^\R) 

















£>[JV-1, 





































V 

where we have also used condition 1 from Sec. 6.2 above in the last step. Now from Eq. (99) 



(E8) 







-)/ L -)/ J 















1^1 J 



(E9) 
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Substituting Eqs. (E8) and (E9) in Eq. (E6) we obtain 

W^, \Wx'f D"' {R) = {R) W^, [W^,f , (ElO) 

where X' is f' or 7', and a and a' are the partitions [N], [N—1, 1] or [A'' — 2, 2]. As we note in Appendix D, 
Eq. (D4), W^,[W^',]^ = for a 7^ a', and so Eq. (ElO) is automatically satisfied when a 7^ a'. It does not 
give us a new condition on the W^/ ■ When a = a' however, we have 

W^,[W^,f D^iR) = D°'{R) W^,[W^,f . (Ell) 

The second part of Schur's Lenima[20] then requires that 

W^,[W^,f = Ia, (E12) 

where is the unit matrix. 

Thus Eq. (E12) is the essential equation to satisfy. 

When Eq. (E12) is satisfied, Eq. (E12), together with the automatically satisfied Eq. (D4), ensures that 
Eq. (E7) is satisfied. Multiplying Eq. (E7) on the right by W'^ and on the left by we obtain 

W^W = I and so W is an orthogonal transformation satisfying 

WW'^ = W'^W = I . (E13) 

F Calculation of [(7^|]x;,x'2 [(7p!^]x'i,x'2 5 the reduced G and FG 
matrix elements, in the symmetry-coordinate basis. 

In this appendix we use the w]^' matrices {X' = v' or 7') to calculate the reduced G and FG matrix 
elements, [cr[^j]x|,x^ and [cr^]x(,x^ , using Eq. (120): 



= (<!)« Qx;x^ [(<')«]" , (Fl) 



where Q = G or FG and ^ is a row label. 
Using Eq. (188) 

{Wf\ = ^\\,.\ (F2) 

and Eq. (197) 



with Eq. (66) 

/ a'liv 
G= I 

g'\M + /I'R'^R 



G=\ . . , (F4) 
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where Ijv , Im , R- are defined in Sect. 5.3 after Eq. (56) and g' and h' are defined by Eq. (65), we find: 

N 

N ^ 

= a', (F5) 



N N k-1 
i=l fe=2 j=l 

= 0, (F6) 



AT j-1 N 

[-[^iW',.' = EEE[<'']».[GT.']..,u(w^i^')"].- 

j=2 i=l fe=l 

= (F7) 

and 

N j-1 N Z-l 

[-r^]]^'.^' = E E E EK"']^. 



j=2 i=i ;=2 fe=i 

W J-1 AT i-1 



2 

EEEE[-'-t"]»j(^''^'^»'=^-?' + ^ii^jk) + h'{6ik + 6ii + Sjk + Sji))^'] 



N(N-l) 

V ' j=2 i=l 1=2 k=l 



kl 



miN -I) 

N 



( E Wi^ik^jl + ^il^jk) + h'{dik + Sil + Sjk + Sjl)] 
i,j,k,l=l 
N 

~ E [d'i^ik^jk + ^ikSjk) + h'{Sik + Sik + Sjk + Sjk)] 

i, j,fc— 1 
JV 

" E [d'i^ik^il + ^iAk) + h'{Sik + Sil + Sik + Sil)] 

i,k,l = l 
N 

+ E [d'i^ik^ik + SikSik) + h'{Sik + Sik + Sik + Sik) 
i,k=l 

[{2g'N'^ - 4h'N^) - 2{2g'N + WN"^) + (2g'N + 4.h'N)] 



2N{N - 1) 

= g' + 2h'{N-l). (F8) 
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Thus we obtain Eq. (201): 



'[N] 



Similarly letting Q = FG and using Eq. (62): 



we can derive: 



FG 



oIn + bJN eR + /J JVM 

cR"^ + d^MN gT-M + /iR'^R + ?Jm 



N 



1 ^ 



ci + bN , 



N N fe-1 
i=l /c=2 j=l 



I — 2 — ^ ^ '^^i 
N V aTTT E E E[1-']' (^('^'^- + + f) tl^'l^-fc 



j=l fe=2 j=l 

=== ( E ^(<5u + Sik) + / - E ^('^'^- + ^^^•) + /) 



= V2(A^-1) (e + y/ 
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N i-1 N 
j=2 i=l fc=l 

1 I — 2 — ^ ^ 

= j^J iy3T E E Ei^^'i'j' r^'^fe' + ^ko) + ^] [if']fe 

* j=2 i=l fe=l 

= ^ /2(N-1) ( ^ ^^^^'^ ^ ^'^'^ + '^l - E [^('^fc* + •^^') + '^l) 



N^2{N-1) 

VW^){<i+^d^ (F13) 



and 

N j-1 N l-l 

K^]h',r = EEEEK^']^.[^^T'T'].i,«[(<'')^].^ 

j=2 i=l (=2 fc=l 



TV i-1 AT i-1 

N{N- 1) 



E E E El-'-'v'liJ [K^ik^ji + '^ii'^jfe) + ^("^ifc + 5ii + 5jk + 5ji) + i\ [l^']fe;) 

j=2 i=l 1=2 k=l 
1 ^ 

■ E {9{SikSjk + Sik5jk)+h{^ik+Sik + Sjk + 6jk)+i) 

i,j,k=l 
N 



2N{N - 1) 



E {di^ikSii + SiiSik) + h{6ik + 5ii + hk + 5ii) + i) 

k,l=l 
N 

E (di^ik^ik + ^ikkk) + h{5ik + Sik + Sik + 5ik) + i) j 

(2iV(7V - l)g + 47V(iV2 - 2N + l)h + N^{N'^ -2N+ 



i,k=l 

1 



g + 2(jV-l)/t+ ^^^^ ^^ i. (F14) 



Thus we obtain Eq. (202): 



[cr^G],,,,-, = (a + bN) = V2(iV-l) + ^/^ 

^ [<^W]h',f^' = VW^) ['^wW,r = (9 + 2{N-l)h+ ^^^- ?j 



(F15) 
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where is the total number of quanta in all the normal modes with the same frequency iv^, i.e. 

oo 

V = d^,n^ . (F17) 



nu=0 
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